Fluctuation-induced noise in out-of-equilibrium disordered superconducting films 
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We study out-of-equilibrium transport in disordered superconductors close to the superconducting 
transition. We consider a thin film connected by resistive interfaces to thermal reservoirs having 
different chemical potentials and temperatures. The nonequilibrium longitudinal current-current 
correlation function is calculated within the nonlinear sigma model description. Different contribu- 
tions are obtained originating from the fluctuation-induced suppression of the quasiparticle density 
of states, Maki- Thompson and Aslamazov-Larkin processes. As a special case of our results, close- 
to-equilibrium we obtain the longitudinal ac conductivity using the fluctuation-dissipation theorem. 

PACS numbers: 74.40.-n, 74.25.F- 
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I. INTRODUCTION 

An equilibrium in nature is rare and it is rather an 
exception than the rule. A vast majority of natural pro- 
cesses, ranging from large-scale flows in the atmospheres 
to the electric charge transfer found throughout the tech- 
nological realm, are out-of-equilibrium processes. Yet, 
the physics of the equilibrium state is much more stud- 
ied and better understood. A major difficulty impeding 
the similar progress in the research of the nonequilib- 
rium processes is that while all thermodynamics-based 
science rests on the law that, in equilibrium, any sys- 
tem assumes the state with the minimal free energy, the 
equally powerful and fundamental principles that gov- 
ern the far-from-equilibrium behaviors still wait to be 
revealed. There has been an important advance both 
theoretical and experimentaliri^ employing a variety of 
approaches to out-of-equilibrium problems in electronic 
systems^"—; among them, those based on the Keldysh 
technique-'^-'iS' seem to appear the most promising as 
paving a way towards general method which would allow 
treating interacting nonequilibrium systems on a com- 
mon fundamental ground. 

In this work we undertake the study of a disordered 
superconducting system employing the theory that has 
been proven to effectively tackle the low-energy excita- 
tion physics, the Keldysh nonlinear sigma model^^. More 
specifically, we focus on the nonequilibrium phase transi- 
tion between the normal and the superconducting state. 
In a vicinity of the transition, when the system is in 
the normal state, the behavior of the system is governed 
by fluctuations of the superconducting order parameter. 
Fluctuation-induced short-living Cooper pairs are formed 
and contribute critically both into the thermodynamic 
and transport characteristics of the systems'^-. In dis- 
ordered thin films the temperature range where fluctua- 
tions are essential is determined by the sheet resistance 
(being in any case significantly larger than the one of 
bulk superconductorsi^"— ) and depends on the particu- 
lar process involved, extending often to temperatures well 



above the superconducting transition temperatureiS."^. 

Here we study the influence of superconducting fluctu- 
ations on dynamic properties of a thin film in the fluctu- 
ational region of a normal state. The film is driven out 
of the equilibrium due to contacts with thermal reser- 
voirs having different temperatures and chemical poten- 
tials. The Keldysh Ginzburg-Landau-like action under 
nonequilibrium conditions and several effects of nonequi- 
librium superconducting fluctuations were addressed in 
Refs.^'^"^^ for this setting. Out-of-equilibrium fluctuation 
contributions to the dc electrical conductivity were calcu- 
lated in Refi^. Importantly, under nonequilibrium con- 
ditions the fluctuation-dissipation theorem (FDT) is gen- 
erally violated and there is no fundamental relation be- 
tween the current fluctuations (i.e. noise) and conductiv- 
ity, in contrast to close-to-equilibrium conditions where 
this relation holds. Therefore, unlike in equilibrium sys- 
tems, the non-equilibrium noise is not fully tied to the 
conductivity and can carry additional information, not 
contained in the conductivity. This poses a problem of 
the independent calculation of noise in out-of-equilibrium 
state. 

Shot noise in noninteracting diffusive electron system 
was studied in Rcfs.^iiJi. The inffuence of the Coulomb 
interaction on shot noise in disordered systems was an- 
alyzed in Refs. Siiiiiii^j while the influence of Bardeen- 
Cooper-Schrieffer (BCS) interaction on shot noise due to 
electric current flow above the critical temperature was 
considered up to the second order in the electric field in 
Refii^. In the present work, we focus on a film above the 
nonequilibrium superconducting transition and calculate 
nonlinear dependence of the Nyquist noise on the tem- 
peratures and chemical potentials of the thermal baths 
that are in contact with the film, see Fig. [TJ 

The paper is organized as follows. In Sec. |ll] we in- 
troduce the nonlinear sigma model for superconductors 
within the Keldysh technique. In Sec. llIII we find the cur- 
rent correlation function for the case of noninteracting 
electrons in a nonequilibrium disordered thin film. We 
further consider the BCS interaction. In Sec.llVlwe intro- 
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duce nonequilibrium fluctuation propagators. Then we 
proceed with calculations of different contributions to the 
current-current correlation function caused by supercon- 
ducting fluctuations: the density of states contribution is 
calculated in Sec.|Vl the Maki- Thompson in Sec. IVII and 
the Aslamazov-Larkin in Sec. IVIII In Sec. IVIIII we sum- 
marize our results. Some calculation details are given in 
Appendices. 



II. KELDYSH FORMALISM: DISORDERED 
SUPERCONDUCTORS 

In this section we introduce the notation and pro- 
vide the basic equations needed for the calculation of 
the current-current correlation function. We present the 
model, discuss its applicability and explain the procedure 
that allow us to analyze the fluctuations of the supercon- 
ducting order parameter in the metallic state. In order to 
treat the nonequilibrium physics, we employ the Keldysh 
techniquei^. We start with the nonlinear sigma model 
and then further discuss the calculation of the current- 
current correlations. 



A. Nonlinear sigma model 

The nonlinear sigma model can be used to describe 
the low-energy physics for superconductors with BCS in- 
teraction in the presence of short-range quenched dis- 
order. The partition function takes the form Z = 
J VQ VAexp{iS[Q]c,^K]}, where the action S consists 
of three part o^^'^^ 



S[Qk,Ak] = Sa + S^ + S, 



Q- 



(1) 



Here and in the following we set h = c = ks — ^■ 
The fields A and Q are Hubbard- Stratonovich fields in- 
troduced to decouple the four-fermion terms originating 
from the Hamiltonian describing the BCS interaction and 
disorder, respectivelyi^^. The contributions to the ac- 
tion IH) are given by 



uses different notation for the same matrices ai = fi for 
convenience, and ctq = diag(l,l). We assume implic- 
itly the multiplication in the time-space, and "Tr" in- 
cludes the integration over the real space. The subscript 
K, denotes the gauge transformed fields 0ac = 4> ~ dtJC, 
Ak; = A + VX; and /C = {k^'ao + k^a^) ® tq. The fields 
A and (j) s^re defined in same way as IC, where A and (f) 
are the vector and the scalar potential, respectively. The 
field A is given by A = {A^'a^ + A^a^) ® t+ - H.c, 
and Ajc(r, t) = exp [ieEIC{r, t)] A exp [— ieS/C(r, t)] . Q/c 
is defined in the same way. We have also defined t± ~ 
{fx ±ify)/2. The quantum (q) and classical (cl) compo- 
nents of the fields are respectively defined as the half-sum 
and the half-difference of the field values at the lower and 
the upper branches of the Keldysh time-contour. The 
field A°' becomes the superconducting order parameter 
at the mean-field (saddle-point) level, while the saddle 
point equation for Q produces the Usadel quasiclassi- 
cal equations, where Q plays the role of the quasiclas- 
sical Greens function. The covariant spatial derivative is 
given by drQic = ^rQic — *e[SAyc, Qk]- We stress that 
the nonlinear sigma model action S captures the low- 
energy physics at energy scales much smaller that the 
elastic scattering rate. It is valid in the limit when the 
lifetime of fluctuating Cooper pairs is much greater than 
the elastic scattering time. 

We further explain the strategy to treat the supercon- 
ducting fluctuations. First we flnd the saddle point equa- 
tion for Q of the action ([1]), in the absence of the BCS 
interaction (i.e. A = 0). It reads asii^i"— 



0=0-2 CiOT^, 

a (r) -U~Hv)-( ~ " 0)^0 Ft,t' (r) 



(4) 



(5) 
(6) 



After Wigner transforming F;/''(r) we obtain F:''\Y,t) 
that can be related to the quasiparticle electron/hole 



-,e/h. 



distribution functions /e^''(r,i) via Ff'^Tjt) = 1 

e/h 



/h. 



5A = -|^Tr[A;^yAd, 



+ iiA,cQK:]- 



Here D is the diffusion coefficient, and it carries informa- 
tion about the disorder. The bare single particle density 
of states at the Fermi level per one spin projection is de- 
noted by h'. The superconductive coupling constant A is 
positive. The matrix field Q satisfies the nonlinear rela- 
tion — 1. The check symbol " denotes 4x4 matrices 
that are defined in the tensor product of the Keldysh and 
Nambu spaces. The former and the latter are spanned 
by the Pauli matrices ai and f^, i G {0,x,y,z}, respec- 
tively, and we define Y = dx ® tq, S = (Tq (E) f^. One 



Tr[0/c^"0/c]i (2) 2/e (r, t). One then considers massless fluctuations 
around the normal-metal saddle point solution, since 
massive modes can be integrated out in the Gaussian 
approximation and lead to unimportant renormalization 
of the parameters in the action. The massless fluctua- 
tions satisfy — 1 and are conveniently parameterized 



(3) 



.17 



:(r) 



A(r) e 



Wir)/2 



WT+ — W*T^ WqTq + WzTz 
WqTq + WzTz WT^ — W*T^ 



(7) 
(8) 



such that WK + KW ~ 0. Here we introduced four 
real fields (r), (r) with a = Q,z representing dif- 
fuson degrees of freedom and the two complex fields 
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FIG. 1. Thin superconducting film is connected by tunneling 
interfaces to tfie gate and tlie substrate. Tlie interfaces are 
cliaracterized by the resistances Ri and J?2. The substrate 
temperature is Ti, the gate temperature is T2 and the gate 
voltage is Va- The parameters Ti, T2, Vg, and the resistances 
of the tunneling interfaces determine the quasiparticle distri- 
bution in the film and allow us to change it in a controlled 
way. 

Wtt'(r), wtt/(r) for Cooperon degrees of freedom. One 
now substitutes Q matrix given by Eq. ([7]) into the action 
(H)) and require that the terms linear in W vanish. This 
leads to a kinetic equation that electron and hole distri- 
bution function have to satisfy. Next we switch on the 
BCS interaction assuming the system is in the normal 
state. In that case the average value of the supercon- 
ducting order parameter is zero, and therefore one again 
obtains the normal-metal saddle point solution (j?]). In or- 
der to study the influence of the BCS interaction, one has 
to consider the fluctuations. Assuming that the system 
is not too close to the transition, we take into account 
quadratic fluctuations around the metallic saddle point 
solution. Integrating them out, one obtains the effective 
action depending only on the superconducting order pa- 
rameter and electromagnetic flelds. This action allows us 
to treat the superconducting fluctuations in the normal 
metalhc state.— 1^ 



B. Current-current correlation function 

Our aim in the following is to calculate the sym- 
metrized two-operator current correlation function 

S{r,t;r',t') = i(j,(r, t)j,(r', t') + j.(r', t')jx(r, t)) (9) 

under nonequilibrium conditions. It can be obtained dif- 
ferentiating the nonlinear sigma model partition function 

b{r,t;r,t) - " 4 5A|(r, i)5A|(r', t') ' a^=o,ac^=o- 

(10) 

In the rest of the paper we focus on the particular sys- 
tem shown in Fig.[T] We consider a thin superconducting 
fllm connected by resistive interfaces to thermal reser- 
voirs having different chemical potentials and tempera- 
tures. We assume that the system is in zero magnetic 



fleld and in the normal state but close to the transition 
into the superconducting state. We study a stationary 

situation (F/^/* = F^/_^,). Therefore, the noise depends 
only on the time difference t^t'. In the following we use 
the notation Ft^^ = F'^/^{e) for the energy dependence 
of the distribution function. 

There are four different contributions to the current 
correlation function 

^(r - r', t - t') ^So{r - r', i - t') + ^oos(r -r',t~ t') 

+ ^MT(r -r\t~ t') + S^^{v -r',t- t'), 

(11) 

where Sq denotes the noise in the noninteracting case 
(A = 0), and the other terms are contributions in- 
duced by superconducting fluctuations. The main pro- 
cesses are: i)coherent Andreev reflections of quasiparti- 
cles on the local fluctuations of the superconducting order 
parameter resulting into the so-called Maki-Thompson 
(MT) contributions-li^; ii) formation of Cooper pairs 
and their involvement into charge transfer is described 
by the Aslamazov-Larkin (AL) corrections^^, and iii) the 
suppression of the single-particle density of states (DOS) 
due to quasi-particle participation in Cooper-pairing—. 
Randomness in the production and dissociation of fluc- 
tuating Cooper pairs, and in other closely related pro- 
cesses discussed above, gives the corresponding contribu- 
tions to the current noise. In the following sections we 
evaluate and analyze the individual contributions. We 
focus on the regime where the superconducting fluctua- 
tions can be treated perturbatively, i.e., the fluctuation 
contribution are small compared to the noise Sq in the 
noninteracting case. This implies that the film is not 
too close to the superconducting transition and that the 
physics is dominated by Gaussian fiuctuations. Then 
one finds S'fluct/5'o ^ Gi = {vDdf)^^ ^ 1. Here Gi 
is the Ginzburg number, rf/ is the film thickness and 
-Sfiuct — S-Dos + Smt + 'S'al denotes the sum of the contri- 
butions induced by fluctuations of the superconducting 
order parameter. 



III. NONINTERACTING ELECTRONS 

We consider noninteracting electrons flrst and evaluate 
Eq. (|10p taking into account nonequilibrium conditions. 
After differentiating the partition function with respect 
to the vector potential fleld in Eq. (|10p , we substitute the 
fleld Q]c with its saddle-point solution A given by Eq. ([4]). 
We do not take here into account the fluctuations around 
A, since they are responsible for weak-localization effects. 
Since the main aim of this paper is to flnd the noise 
close to the superconducting transition, and since the 
weak-localization contribution is a non-singular function 
in the vicinity of the transition, we do not consider it 
here. Then, we obtain the longitudinal current-current 
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correlation function to beii 

So{r~v',e)^2<jDTo{e)S{r~v'), (12) 

+ F\n)F\n + e)]y (13) 

Here the Drude conductivity is aj) = 2e^vD and the 
integration is always from — cxd to +oo if not specified 
differently. In equilibrium, the distribution function is 
F'=/''(e) = tanh(e/(2T)) and the expression ()12p simpli- 
fies to 

Sl\v - r', e) = 5{v r')az,ecoth (^) . (14) 

This universal relation between the conductivity and the 
noise is known as the FDT. The microscopic details about 
the disorder strength are hidden in the diffusion constant, 
i.e., in an- We stress that the range of applicability of 
the nonlinear sigma model is e ^ t~^, where r is the 
elastic scattering time. This explains why in Eq. (fT4|) 
appears the frequency independent Drude conductivity 

Next we analyze Eq. (|12l) in out-of-equilibrium con- 
ditions. Then, the FDT is generally violated. We as- 
sume that the film is thin such that the Thouless energy 
corresponding to diffusion across the film E;^ = D/dj, 
well exceeds all the relevant energy scales. Here df de- 
notes the film thickness. The current across the in- 
terface separating the substrate and the film is / = 
/ de [F^ie] - T|(e) - F''(e) + F^{e)] /(4ei?i), where Ri 
is the tunneling resistance per unit area characterizing 
the interface and the subscript S denotes the substrate. 
A similar equation holds for the interface between the 
film and the gate. From the continuity equation for 
the current follows F'^/''(e) = xFg^'\e) + {1 ~ x)FG'^{e), 
where x = i?2/(i?i + i?2)- Here Fl'^{t) = tanh(e/2ri) 
and Fa'^{e) = tanh ((e =F eVG)/2T2) denote the distri- 
butions in the substrate and in the gate, respectively. 
The gauge invariant distribution in the film is defined as 
F'=/''(e) = F'=/''(e ± e(tf^) and takes the form 



xtanh 



e±{l-x)eVG 



+ {1 — x) tanh 



2Ti 

e =F xeVc 
2To 



(15) 



The upper (lower) signs correspond to electrons (holes). 
We assumed very resistive interfaces, such that the re- 
sistance of the film can be neglected with respect to the 
resistance of the interfaces. 

Now we can proceed with the evaluation of the expres- 
sion (fT2|) . In the case Ti — T2 — T one can calculate it 
exactly and obtains 

Soir - r', e) ^S^^v - r', e)[x' + (1 - xf] 

+ S^''{r^r\e + eVG)xil~x) 

+ S^\r-r',e-eVG)x{l-x), (16) 



where Sq'' is the equilibrium noise, see Eq. (IT4l) . When 
one of the tunneling resistances is infinitely large, then 
effectively the system is in contact only with one reservoir 
and therefore in an equilibrium. In that case a; = or 
X = 1, and as expected Eq. (fT6)) reduces to the equi- 
librium noise Sq"^ . Moreover, notice that by increas- 
ing the temperature of the system, thermal fluctuations 
increase and the current noise increases. Similarly, in 
Eq. ((T6)) . the gate voltage also leads to an increase of 
the noise with respect to the equilibrium one {Vg — 0). 
In the case of zero frequency, the noise 5*0 (r — r', 0) was 
found in Ref.^**. Eq. pB]) also describes a diffusive bridge 
placed between two reservoirs having different chemical 
potentials'^. There, x plays the role of the coordinate 
along the bridge, and after performing the integration 
over X in Eq. one finds the shot noise of that sys- 
tem. 

Next we discuss the influence of the electron-electron 
interaction on the distribution function (|15p. In general, 
the interaction leads to smearing of a noninteracting dis- 
tribution function. However, when the inelastic length 
is large in comparison to the system dimensions, inelas- 
tic processes can be neglected. The inelastic length is 
expected to increase as the temperature and/or applied 
voltage decrease, allowing one to tune the ratio of the 
characteristic system length and the inelastic relaxation 
length. Coulomb-interaction induced corrections to the 
noise in the equilibrium have been studied by Altshuler 
and Aronov"^^. Out-of-equilibrium, the influence of the 
the Coulomb interaction on the shot noise in diffusive 
contacts has been studied in the limit of large inelastic 
length in Refs3 '^^i^^ , while contributions due to inelas- 
tic collisions have been studied in Refi^. In the present 
paper, we assume that time-of-flight of the quasiparti- 
cles through the film is much smaller than the typical 
energy-relaxation time. We consider a thin film charac- 
terized by the large inelastic scattering length and focus 
on the influence of the BCS interaction on the longitudi- 
nal transport in the normal state. 



IV. SUPERCONDUCTING FLUCTUATIONS 

Having discussed the noninteracting case, we start the 
analysis of the influence of the BCS interaction on the 
current noise. The system is assumed to be in the nor- 
mal state, but in the vicinity of the nonequilibrium su- 
perconducting transition. The saddle point equation of 
the action ([1]) for the superconducting order parameter 
has the trivial solution with the average value of the su- 
perconducting order parameter being zero. Then, one 
recovers the normal-metal saddle point ^ for the field 
Q. Therefore, now we include the massless fluctuations 
around it. We substitute Eq. ([7]) into the action (P) and 
expand it to the second order in W, Eq. ([8]). Since we 
are neither interested in the weak localization correc- 
tion nor in the Altshuler-Aronov type corrections, but 
in fluctuations caused by the fluctuating superconduct- 
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ing order parameter A(r,t), in the following we consider 
only the Cooperon degrees of freedom. After integrat- 
ing out the Cooperon degrees of freedom, we obtain a 
Ginzburg-Landau-like action which depends only on the 



order parameter A{r,t). This calculation was performed 
in KeiM^ and the correlation functions of Cooperon fields 
have been obtained. They read as^,— 



/ \ * r \\\ 2i —LJ^LA,i-2LR^i-2 + F'^{e3)Ljii-2 + F'^{^i)LA,i~^2 



(Kl,e2(q)w^-e3,-e4 

(('«ei,e2(q)w^-e3,-e4 



2i 

(q))) = —Sei-e2,e^-es 



— Lj}^La,1-2Lr^I-2 — F^{(-2)La,1-2 " -F'^(e4)ifl,i-2 
[i?q2 + z(ei + 62)] [Dci^ + iie-s + 64)] ' 

Lj(^La,1-2Lra-2 + F^{^2)La,1-2 — F'^{^?,)L]i,i^2 

W+i(ei+e2)] W- 1(63 + 64)] ■ 



(17) 
(18) 
(19) 



The average ((. . .)) is with respect to the action S given 
by Eq. ([T]) and it includes averaging over the fluctuations 

of Q, A^' and A«. Here Lr/a.i-j = (-^fl)yi(q: e^ - Cj)) 
denotes retarded/advanced fluctuation propagators and 
is the Keldysh propagator. Their low frequency w 
and low momentum q behavior is given by the following 



expression! 



.25. 



^ 2 



1 -F''(0)F'=(0) 

TT f 



(20) 



+ 



T 4iT^F'^{0) - Dq^ ± ie T iieVail - x) 
T 

(21) 



(-'If)}' 



Here and in Eqs. (I17m9|) we set Ajq = 0. The parame- 
ters appearing in the retarded and advanced propaga- 
tors are functional of F"^{e) = [F''(6) - F'=(-6)] /2. 
They are given by the following expressions T^^ = 
2d,F^(e) , f^-i = 2f d6 [/'«(6) -F^(0)1 /(627r), 

and = 1-1- (Te/ri) . The symbol / denotes 

the principal value of the integral. The nonequi- 
librium Ginzburg-Landau (GL) rate is defined as 



de 



tanh(6/2Tc) /ire + 

^F^{0), where ujd is the Debye energy. The 
GL time denotes the lifetime of the fluctuation induced 
Cooper pairs. It carries the information how far the sys- 
tem is from the transition and at the transition it become 
infinitely large, signaling that the Cooper pairs become 
long living. We point out that the fluctuation propa- 
gators given above are valid when the system is in nor- 
mal state but in the vicinity of the transition such that 
(tglT^c)"^ ^ 1- Substituting the distribution function 
(1151) in the previous equations, one evaluates all these 
parameters and finds that they are functions of temper- 
atures of the reservoirs Ti and T2, the ratio of tunneling 



resistances x — i?2/(-Ri + R2) and the gate voltage Vq^: 

r 1^^ 

J. _ X ^ {l-x) 



Ti cosh^ T2 cosh^ ^ 



2x 



rim 



2{l-x) 



^, (l_^ eVM-x) 
2 2ttTi 



Im 



2 27rT2 



(22) 



(23) 



-ZcpTe< xRe 



(1 - a;)Re 



^ I il-xWo 
2 2ttTi 



^ , I , .xeVc 
2 27rT2 



21n2 



+ +4zcp§F^(0). (24) 



Here "^(z) is the digamma function, defined as 
^(x) — r'(x)/r(a;), where r(a;) is the gamma func- 
tion. We expressed the critical temperature as Tc = 
2a;£) exp (7 — A~^)/7r, where 7 is the Euler constant. Now 
one easily finds the parameter z~p =1-1- (Te/Sl)^. 

We are now equipped to start the calculation of differ- 
ent fluctuation contributions to the current-current cor- 
relation function. However, before doing it, we make 
a brief digression in order to stress the importance of 
the above given fluctuation propagators. In close-to- 
equilibrium conditions, the Aslamazov-Larkin contribu- 
tion to the conductivity can be obtained by employing 
together the linear response theory and the phenomeno- 
logical time-dependent GL (TDGL) equation, while the 
Maki- Thompson and the DOS, should be derived start- 
ing from the microscopic theory. The phenomenological 
TDGL equation reads as 



(T - T,) - D{ci - 2ek^f + i{uj - 260^) 



A^(q,c.) 



+C(q,c^) = 0, 

(25) 
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where the thermal noise satisfies 
1 fi 

{C{r,t)C{r',t'))^—T'S{r~r')5{t~t'). 
in/ 



(26) 



However, if one is interested in a nonhnear dependence 
on the drive, then the previous equation is not a good 
starting point. Notice that even for the simple setup 
here considered (see Fig. [T]), the TDGL equation takes 
different form^^ than the usual phenomenological TDGL 
equation ([25)) . It is given by 



L^i(q,w)A^(q,c.) + C(q,^)=0, 



(27) 



where the nonequilibrium noise satisfies the following 
condition 



(C(r,t)r(r',t')) 



16, 



1 -F''(0)F^(0) 



X 5{T-Y')5{t-t'). 



(28) 



In Eq. (f27| . the retarded fluctuation propagator is given 
by Eq. pi|) with many drive-dependent parameters that 
we discussed above. Only in the equilibrium, Vg = 
and Ti = T2, Eqs. ((271) and ^ coincide with Eqs. ^ 



and (|26p. while otherwise they are different. Therefore, 
we point out that in order to study a nonlinear depen- 
dence of some observable on the drive, one should derive 
a corresponding TDGL equation starting from the micro- 
scopic theory and not to use the phenomenological one 
([25| . as is usually the case in literature. 



V. DENSITY OF STATES CONTRIBUTION 

In this section we start the calculation of the current 
correlation function ([T0|) in the presence of the BCS inter- 
action. We consider linear response to the in-plane elec- 
tric field, while the system is driven out-of-equilibrium 
due to contacts with thermal baths at different tempera- 
tures between which it is sandwiched, or due to an electric 
field perpendicular to the plane, see Fig. [TJ 

First we consider the density of state contribution to 
the noise. Taking into account the massless fluctuations 
around the normal-metal saddle point Eq. (j4]) to the sec- 
ond order in the Cooperon fields, we collect all the terms 
of the type {{ww*)) and {{ww*)) in Eq. dTU]). They con- 
stitute the DOS contribution and together give 



S'oos(i" - r',e) = 



~ Wit 

X [1 - F'^ie^ ~ e)F\e2) - F^{t^ - e)F^{e,)\ - {{w,,,,A^)w 



5{T~r') / deide2<j ((Wei,e2(i-)w*e2,-ei(r) +'«n,f2(r)^«* £2,-61 W)) 



xF''{e2)F^{ei+e)^+e^-e. 
The leading contribution in Eq. (1291) close to the transition into superconducting state is given by 



(29) 



5oos(r-r',e)= + (5(r-r')Y^^ j d\AEAujlui(^-^^ La^uLb..u^2 



F'^{E - e)F^{E) 
' [Dq^ ~ i{2E + uj)Y 



F^{E -uj)F''{E + e) 



[Dq^ - i{2E + 2e - w)] [Dq"^ - i{2E - to) 



(30) 



We use the notation L 



R/A,i 



i^J^(q, w) ) . Notice that Eq. (UHl) is the only singular part of Eq. for r^^ ~^ 



and therefore the most dominant close to the transition. The holes and electrons give the same contribution. This is 
expected, since the system is in the normal state and there is no long-living Cooper-pair condensate that could allow 
for a charge imbalance between the holes and electrons. We further evaluate Eq. ([501 and find up to logarithmic 
accuracy: 



S'oos(r - r',e) 



27r2d 



-5{v-y') 



1 -f'''(0)F^(0) 



Zcp In 



Re< 2 / dE 



^F'^{E - e)F'^{E) 
{E + i0)2 



dE 



F^{E)F%E + e) 
{E + € + iO){E + iO) 



e -!• -e. 



(31) 



Here and in the following iO denotes ?0^. In order to discuss the assumptions made in Eq. pip , we introduce 

eo = -4TeF«(0) + T,T-^/n. (32) 
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We assume that the system is close to the transition such 
that I eg I and are much smaller than relevant energy 
scales of the distribution function, i.e., min(VG, Ti, T2). 
For the special case x — 1/2 and Ti = T2, eg is exactly 
zero. Eq. ((3T|) is also valid for Vg = and then the con- 



dition becomes |eo|,TQL <C min(Ti,T2). When obtaining 
Eq. pil) . we used the fact that the most important con- 
tribution in Eq. (j30p originates from small momenta and 
that therefore we can safely cut the momentum integra- 
tion at the upper limit Dq^^^ ~ Tg. 

In the case of equal temperatures of the reservoirs Ti — 
T2 = T, substituting the distribution function ([15]) in 
Eq. (|3T|) . we obtain the nonlinear dependence of the noise 
on the gate voltage Vg, temperature T and frequency e. 
Integrating out spatial coordinates, one obtains the final 



result 



'S'dos (e) 



1 -/'''(0)F"(0) 



TeZcp In ; 

Tr.T, . 



x^D{eVG{l-x),eVGil-x)) 

+ x{l - x)D (eVfe(l - x), -eVcx) 
+ x{l - x)D (-eVcx, eVcil - x)) 

+ (1 - x)^D {-eVcx, -cVgx) + e ^ -e 



(33) 



The function D(x, y) is defined as 



D{x,y) =-coth 



y — X + e 



2T 



Im 



^ \ I 

2 27rr 



y 



2 2ttT 



y 



7^°*M-2T 



X Re 



2ttT 



X 



2ttT 



1 . y 

h t 

2 27rT 



■ y 



2ttT 



(34) 



where the frequency satisfies e ^ |eo|,T(7L. However, 
for e <C min(VG, Ti, r2) but arbitrary with respect to 
|eo|,T~L , one obtains the result by setting e to zero in 
Eq. p3p . The expressions for Tql, Tg, Zcp and ft are 
provided in Sec. IIVI In Appendix |^ we discuss the case 
of different temperatures of the reservoirs, Ti ^ T2. 

Next we analyze the expression (1551) . For a; = (or 
a; = 1) the system is effectively in an equilibrium and 
Eq. ([55)) simplifies to 



2fT 

TT^df 



■In 



T 



coth ( 

V2r 



Im 



--i — 
2 2ttT 



T 

27r-Re 
e 



- ^ 



2 ^2ttT 



(35) 



For Vg = we also get the same expression from Eq. 
since it also corresponds to the equilibrium situation. In 
the case of noninter acting equilibrium electrons, stud- 
ied in Sec. IIIH the source of the current fluctuations is 
thermal fluctuations. Here we obtained the additional 
contribution ()33p to the current noise caused by super- 
conducting fluctuations. The Cooper pair density fluctu- 
ates due to randomness in the formation and dissociation 
of Cooper pairs, and affects the single-particle density of 
states leading to the noise (S'dos- 

In equilibrium the FDT holds and relates the current 
fluctuations (noise) to the real part of the ac conductivity 
that provides information about absorbed energy in the 



sample: 



S"=«(e) 



Re[a(e)]ecoth(^), 



(36) 



where S{e) — J dxS'(x, e). Using it we obtain the density 
of states contribution to the in-plane ac conductivity in 
equilibrium: 



2e^ T 
TT^df e 

T 

27r-Re 
e 



21e"C(3) 



In 









m 


Im 











2 *27rT 



2ttT 



(37) 



In 



-J^ln{^){^)'ln{^), e»T. 



(38) 



This result is in the agreement with findings of Refj^, 
where it is obtained in a different manner. In 
Eqs. P5I37|38|) the GL rate takes its equilibrium value 
T~j^ — STln (r/Tc)/7r. One obtains that suppression of 
the single-particle density of states due to superconduct- 
ing fluctuations gives the negative contribution to the 
Drude conductivity, see Fig. [2] where S'dos as the func- 
tion of frequency e is shown by the red dashed line. We 
remind the reader that our calculations apply when the 
frequency, temperatures and the gate voltage are much 
smaller than the elastic scattering rate, since we use the 
nonlinear sigma model that captures low-energy physics. 
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5 10 15 

FIG. 2. Current noise Suos{^) is shown as a function of fre- 
quency e for the case of equal temperatures of the reservoirs 
Ti — T2 — l.OlTc, the ratio of tunneling resistances x — 1/3. 
The red dashed line presents zero gate voltage and the blue 
solid line Va = Tc/2. 

Now we consider the influence of the gate voltage on 
the noise S'dos- From Eq. ([55)1 one easily finds that the 
first correction to the equilibrium noise (j35p in the limit 
Vg <^T depends quadratically on the gate voltage. This 
is expected since in the case of equal temperatures of 
the reservoirs, Ti = T2 = T, in-plane current is invari- 



ant under the transformation Vq —Vg and therefore 
S'dos, Eq. ([55)1 . is the even function of the gate voltage. 
Fig. [2] shows S'dos, given by Eq. (|33l) . as a function of e 
for the case Ti = T2 = 1.0 ITc and the ratio of tunneling 
resistances x = 1/3. The red dashed line represents zero 
gate voltage, i.e. the equilibrium case, while the solid blue 
line is for Vg = Tc/2. One observes that the gate voltage 
tends to suppress (the absolute value of) the equilibrium 
noise, Sdos, contrary to the noninteracting electron case 
where it increases the noise, see Sec. IIIII The reason is 
that the gate voltage decreases the Cooper pair lifetime 
and their density, and therefore the infiuence of supercon- 
ducting fiuctuations on the current correlation function 
decreases. The lifetime is given by the GL time Tgl, see 
Eq. 



VI. MAKI-THOMPSON CONTRIBUTION 

In this section we consider the Maki-Thompson con- 
tribution to the current-current correlation function. Its 
physical origin is already discussed in Sec. |ll] and here 
we calculate it. We start from Eq. and collecting 
all the terms that contain the average of the convolu- 
tion of w and w* fields, we find the the Maki-Thompson 
contribution 



SMT(r - r', e) ^^S{r - r') J dEduj{{{wE+.Mr)wlE,,^,E_,^Jr))) [-1 + F\E)F\E + e)] 
+ {{wE+^A^)w*_E_^^_E_^_^{v))) [-1 + F'^iE + ^)F'^{E + e + 0.)] + e ^ -e}. 
The leading contribution of the previous expression close to the superconducting transition reads as 



5*mt(^) 



2-Kdf 



I- F^{Q)F''{Q) 



2r, 



{2 -F''(0)[F''(e) 



(39) 



(40) 



and is valid for i) e < min{Ti, T2, Vg} or ii) Vg = 
and e ^ minjTi, 12}, but arbitrary with respect to r^^ . 
The condition for Tgl and eo, Eq. in comparison 

with the temperatures and the gate voltage is the same 
as for the density of states correction, see the discussion 
below Eq. (|32p . For the definitions of the parameters 
appearing in Eq. (HUl) . see Sec. IIVI while the distribution 
function is given by Eq. (jlSp . One sees that S'mt is not 
singular function of the GL rate at finite frequency e. 
In the case e ^ r^^ ^ t"^ , the phase breaking rate 

appears in Eq. PO)) inside the logarithm instead of 
e, because it serves as a cutoff scale at small momentum: 

Note that contrary to the density of states contribu- 
tion where the characteristic frequencies are determined 



by the temperature and the gate voltage, here two addi- 
tional energy scales appear: the GL rate and the phase- 
breaking rate. Fig. [3] shows Smt as a function of fre- 
quency e for the case Ti — T2 — l.OlTc and the ratio 
of tunneling resistances x = 1/3. It is assumed that e 
is greater than the phase-breaking rate. The red dashed 
line presents equilibrium case (zero gate voltage) and the 
blue sohd line corresponds to Vg = Tc/2. We see that 
Smt monotonically decreases with frequency, and is pos- 
itive contrary to the negative Sdos contribution. On the 
other hand, the role of the gate voltage is similar in both 
Sdos and Smt, it Suppresses superconducting fiuctuations 
and therefore their infiuence on the current noise. 

Further, we discuss the infiuence of temperatures of 
the thermal baths on the noise. For simplicity we con- 
sider the zero gate voltage. The system is driven out- 
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FIG. 3. Current noise 5'MT(f) cis 3- function of the frequency 
e ^ for the case of equal temperatures of the reservoirs 
Ti = T2 — l.OlTc and the ratio of tunneling resistances x = 
1/3. The red dashed line presents zero gate voltage and the 
blue soUd Une Va = TJ2. 
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ferent realizations. The ratio of tunneling resistances is 
fixed to be x = 1/3. The red dotted curve corresponds 
to the equilibrium situation Ti = T2 = l.OSTc- Lower- 
ing one of the temperatures, the effective temperature 
of the system decreases, the lifetime of the Cooper pairs 
increases and therefore the noise also increases, as can 
be seen by comparing the blue solid curve in Fig. 21 that 
corresponds to the experimental realization Ti = O.OSTc 
and T2 — l.OSTc, with the red dotted curve. On the other 
hand, by increasing one of the temperatures, the effec- 
tive temperature of the system increases, the lifetime of 
Cooper pairs decreases and therefore the influence of the 
superconducting fluctuations decreases, leading to the re- 
duction of the noise, as shown by the black dashed curve 
for the case Ti = 1.15Tc and T2 = l.OSTc. 

In equilibrium, Eq. (j40[) further simplifies. Using the 
FDT we find the leading MT ac conductivity^^ 



Re[crMT(e)] 



■2t, 



2nd 



-T- 



/ 

■_ef_ T 
T 



To.' 



(41) 



ln(r-,7|e|), T » r » e 



T > e > r: 



(42) 



Here the GL rate takes its equilibrium value 
8rin(T/Tc)/7r and r^^ < e < T. 



In the case e 
-1 



^ ''gTl I the phase breaking rate appears in in- 
side the logarithm instead of e, because it serves as a 
cutoff scale at small momentum: Dq^^-^^^ ^ '''s^ ■ 



VII. ASLAMAZOV-LARKIN CONTRIBUTION 



FIG. 4. Current-current correlation function Smt as a func- 
tion of frequency e for different temperature of the baths. The 
ratio of tunneling resistances is fixed to be a; = 1/3 and the 
gate voltage is Vg ~ 0. 



of-equilibrium due to different temperatures of the reser- 
voirs. Fig. 2] shows frequency dependence of S'mt for dif- 



In this section we obtain and analyze the Aslamazov- 
Larkin contribution to the current-current correlation 
function. Having discussed the density of states and 
the Maki-Thompson terms in Eq. (fTOt . all the remain- 
ing terms are the fourth order in the Cooperon degrees 
of freedom. They constitute the Aslamazov-Larkin con- 
tribution. We use the notation S'AL(e) = / drS'AL(r, e) 
and find 



5'AL(e) 



1 /nveDV 



ei,e2,e4,e5,qi,q3 



F''(ei)u;,„,,_,(qi)iDl,^,_,^(qi) 



+ i^''(ei)wl,,,_,,+e(-qi)wei,e2(-qi) + ^''(ei - e)we,,e,-,{qi)wl,,^-,,{qi) 

- F'^iei - e)w;*,^ _,^+,(-qi)w;,,,,,(-qi) - F''{e4)we,,e^+eiq3)wl,^,^eM3) 
+ F'=(e4)wl,^^_,^_,(-q3)we4,c5(-q3) + P'^ie^ + e)w,5,,^+,(q3)wl,^__,^ (qa) 

- F^{e4, + e)w* ,,,_,,_,(-q3)we4,e5(-q3) 



(43) 
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where = J dei/(27r) and J = J dqi/(27r)^. The leading contribution in Eq. p3)) close to the transition into the 
superconducting state is given by 



1 r meDT^Zcp 



Im 



dE- 



1 -F''(0)f"=(0) 



2 /.+C 



dco 



Vg ^ ~Vg 



+ 00 



dqq^ 



(tgl + Dq^)^ + (e - 0^)2 (tc,^ + Dq^)^ + lo^ 



(44) 



The condition of applicability of the last formula is the same as for the density of states correction, see the discussion 
below Eq. ([32|) . Since we used the fluctuational propagators from Sec. |lV]that are applicable for low frequencies, 
Eq. (j44l) is valid for e much smaller than the reservoir temperatures and the gate voltage. However, Eq. (|44p is valid 
also for Vg = and then the condition is e ^ min(Ti,T2). Evaluating Eq. (j44|) . we find 



l-F 



''(o)i^"(o)]y 



■ arctan 



■ In 1 



-Re 



27rTi 



il-x)eVG 





l-x 




)] 








^2 





27rT; 



4ro 



■xcVg 



Vg ^ -Ife 



(45) 
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FIG. 5. Current noise Sai,{^) as a function of frequency e for 
the case of equal temperatures of the reservoirs T\ = T2 = 
l.OlTc and the ratio of tunneling resistances x — 1/3. The 
red dashed line presents zero gate voltage and the blue solid 
line Vg = n/2. 
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FIG. 6. Current-current correlation function Sa-l as a function 
of frequency e for different temperature of the baths. The 
ratio of tunneling resistances is fixed to a; = 1/3 and the gate 
voltage is Vg = 0. 



T~} and Te see SecHVl The 



For the definitions of 
distribution of electrons and holes in the film is given by 
Eq. (IT5|) . where x = i?2/(^i + ^2) is ratio of tunneling 
resistances, and d/ is the film thickness. 

Note that the only relevant energy scale of the fre- 
quency dependence of Sal is the GL relaxation rate. 
Fig. [S] shows Sal as a function of the frequency e for 
the case Ti = T2 = l.OlTc and the ratio of tunnel- 
ing resistances x — 1/3. The red dashed lines presents 
zero gate voltage and the blue solid lines corresponds to 
Vg — Tc/2. One sees that the Aslamazov-Larkin contri- 
bution is strongly suppressed by the gate voltage due to 
suppression of the Cooper pair lifetime. The influence of 
temperatures of the thermal baths on the noise is shown 
in Fig. ini One observes a similar behavior as in the case 



of the Maki Thompson contribution. 

Finally, in the equilibrium (Ti = T2 ^ T and Vg = 0) 
we find the leading term in the AL ac conductivity using 
the FDT and Eq. ^^^^ 



Re [cT Ai^ie)] =— 



2f T 
ndf e 



arctan 



GL 

e 



In 1 



GL 



4tgl' 



(46) 



Here the GL rate takes its equilibrium value t^^^ — 
STln (T/Tc)/7r. We point out that subleading terms to 
Eq. in the dc case logarithmically depend on the GL 
rate, see Appendix IB] or Ref."^^ where they are obtained 
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in a different way. 
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Appendix A: Density of states contribution to the 
noise 

In this section we discuss the density of states contri- 
bution to the noise for the case Ti ^ T2, close to the 
transition into the superconducting state. In Sec. |V] we 
calculated it for Ti = T2, see Eq. ([55)) . We start form 
Eq. ([3T|) and obtain 



'S'dos(£) 



n^df 



TeZcp In ( ^ 

'''gl 



x^D, ieVG{l-x),eVGil-x),T,,T,) 

+ x{l-x)Di {eVG{l-x),-eVGX,Ti,T2) 
+ x{l-x)Di {-eVGX,eVG{l-x),T2,Ti) 
+ (1 - xfDi i-eVGX, -eVGX, T2, Ta) 

(Al) 



Here the function Di is defined as 



Di{x,y,Ti,T2) =4Im 



1 tanh(l^+^2l^7rf 



tanh 



Ti , x+y+€ 
T2 2T2 



1 ^tanh(m+,2^,|^ 



2Ti 

-E 



^ tanh iTT^f -^if±i 



n>0 



^(2n+l)-^ 



Tl>0 
1 



(2n+l + z^)2 



z(2n+l) + 2i^ z(2n + l) + -f^ 



z(2n+l)-|^j 



(A2) 



For the validity of this expression see the discussion below Eq. ([5^ . The definitions of the parameters are given in 
Sec. IIVI From Eq. (jA2p . we obtain the closed-form result for large frequencies, i.e. in the limit |e|, |e ± eVcj ^ Ti, 12 



Di{x,y,Ti,T2) =sgn(2/ - a; -1- e)Im 



1 T/ / l .~x + e\ 1 /I . y 

— t * - + i—^ 

Ti I 2 2ttTi I T2 I 2 27rT2 



— sgn(2: + y + e)Re 



T I 1 ■ ^ 

I 

2 2ttTi 



2 27rT2 



2 27rT2 



T I 1 .x + e 

* I 

' 2 2ttTi 



Q(^Q-\y-x+<^\/Ti^^-\v+x+e\/Ti^^ 



(A3) 



where i = 1,2. In the case Ti ^ T2 = T, the func- Eq. ((M)) . 
tion Di{x,y,T,T) becomes equal to D{x,y) defined by 
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Appendix B: AL conductivity and using the FDT. We obtain the leading term close to 

transition to the superconducting state to be 

In this section we calculate the Aslamazov-Larkin ac 
conductivity in the equilibrium, starting from Eq. (j43p 



oil) e^D 



2 n2 r+oo 



87r2 

4-CXD 



/ dqq' du\L^{q,io)f\L''{ci,LU-e)\'L],\q,u)L],\q,u~e) 

Jo J ~oo 

tanh 



dei 



tanh (l^j 



[L»g2 _,_ i(^2ei + w - 2e)] [Dg2 + j(2ei + a; - e)] [1592 + ^(2^^ _ a;)] [Dg2 + i(2e^ - w - e)] 



(Bl) 



In the following we use the complete expressions for the fluctuation propagators, since they are needed in order to 
calculate the subleading terms. The fluctuation propagators read as"^*^ 



1 _Dq2 — iuj 



AttT 



-^K^(q,'^) = coth (^^) [L^/(q,w) - L^\q,uj)] 



where the Ginzburg-Landau rate is defined as 



(B2) 

(B3) 
{B4) 



-1 ST T, 
^GL = In I — 

TT 



T 



(B5) 



Using 



+°° , tanhx 2 
dx = 

,00 {x + a){x + b) a — b 



1 ia 



1 ib 



Im(a),Im(6) < 0, 



(B6) 
(B7) 



we perform the integration over ei. Then, analyzing this expression, we find that the main contribution comes from 
small momenta and for frequency |a;| ^ and |w — e| ^ . In the limit t^T^ , e <g T but for an arbitrary ratio of e 



and Tqj^, we find 



;.(1,2) 



Tidf e 



2 GL 



In 1 



14e2C(3) T 



4eTrr, arctan 



4 

+ ln( 1 



TglC arctan 



£2x2 



2 2 1 



4 + eV2, 



(B8) 
(B9) 



Here the first line gives the leading and the second one the 
subleading contribution, S'ii''(e) ~ S^Ajl^\e) + 5'iL^^(e). 
Eq. (j45p coincides with sIl^'' in an equilibrium. 

We calculate further an additional contribution to the 
AL noise, originating from Eq. (|43p. and not taken into 
account in Eq. (jBll) . It is also singular function for 
T(7l ^ in the dc limit, but less relevant than the ex- 
pression (jB8l) . It gives a contribution to the noise that 
is of the same order as the subleading term (IB9p . The 
Aslamazov-Larkin noise, Eq. (j43p . contains products of 
two correlation functions given by Eqs. (|17m9l) . In the 



leading contribution (jBll) . when calculating correlations 
(I17m9p . only the terms ~ were taken into account. 
In the present case, when calculating product of the cor- 
relation functions, we should take into account the prod- 
uct of one term ^ and one that contains only a 
retarded or advanced fiuctuation propagator. Then we 
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find the leading term 



n^df 



■Im(/i)Re(/2) 



+ 00 



dqq^ 



+ 00 



duj\LA{q,uj)\'^Lj/{q,uj) 



xRe[LA{q,Lo-e)+LAiq,Lo + e)]. (BIO) 



tions Ii and I2 are defined as 

^+00 



Ii=- 



tanh(e/2T) _ . tt 
+°° tanh^(e/2T) _ 7C(3) 



Evaluating the expression (jB10[) we obtain 



c(2,l) 



112C(3)e2^T^ 



4Tn 



■In 



arctan 



2r„ 



y2 



+ £2/4 



(Bll) 
(B12) 



(B13) 



Here the lowest order term in the expansion of the fiuctu- 
ation propagators for small uj and q appears. The func- 



Therefore, in the limit t^l ,e <^ T we obtain the lead- 
ing term S'li'^'' given by Eq. (jB8l) and the next leading 

fl 2^ ('2. 1) 

order is the sum of S'Il and Sal given by Eqs. (|B9|) 
and (|B13|) . respectively. Then, using the FDT, we find 
the first two leading terms in the real part of the AL 
conductivity 



Re [ctal (e)] 



2e^T 
^-1 



arctan 



2Tn 



^In 

4e 



4Tn 



-i^ln 
e 



4to 



■In 



4roL' 

r2 



28C(3)e 



arctan 



TGif + eV4 



2Tn 



(B14) 



The leading term is in the agreement with the calcula- 
tion of Refsj^i^^, and the subleading one was recently 
obtained in Reffi'* in a different way than the one pre- 
sented here. From Eq. (IB14[) . in the dc case we obtain 

r . M T 14C(3)e2 , / T \ , 

ReK.O =— -— + ^Ij In — . B15 



Apart from the well known power law singularity of the 
Aslamazov-Larkin dc contribution for t^-^ — > 0, we get 
the additional logarithmic one in agreement with Refi^. 
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